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We show that there is a class of stochastic baker's transformations that is equivalent to the class 
of equilibrium solutions of two-dimensional spin systems with finite interaction. The construction 
is such that the equilibrium distribution of the spin lattice is identical to the invariant measure 
in the corresponding baker's transformation. We also find that the entropy of the spin system is 
up to a constant equal to the rate of entropy production in the corresponding stochastic baker's 
transformation. We illustrate the equivalence by deriving two stochastic baker's maps representing 
the Ising model at a temperature above and below the critical temperature, respectively. We cal- 
culate the invariant measure of the stochastic baker's transformation numerically. The equivalence 
is demonstrated by finding that the free energy in the baker system is in agreement with analytic 
results of the two-dimensional Ising model. 
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The Ising model and the baker's map [l| are two exam- 
ples of canonical model systems that have provided deep 
insights into statistical mechanics and dynamical systems 
theory, respectively. The aim of this paper is to point 
out and illustrate a relatively simple connection between 
spin systems in two-dimensions and a class of stochastic 
baker's transformations. In this way we demonstrate a 
new connection between equilibrium lattice systems and 
low-dimensional dynamical systems. 

Formal connections between dynamical systems and 
two-dimensional spin systems have been presented be- 
fore. Examples include cellular automata (e.g., 0) and 
coupled map lattices (e.g., [1, 13]). We introduce a new 
class of stochastic baker's transformations and show that 
the noise can be constructed so that the representation of 
the dynamics is identical to a representation of equilib- 
rium spin systems in two dimensions Q . The novelty in 
the present approach is the design of a specific family of 
transformations in which the invariant measures are iden- 
tical to the equilibrium distributions of two-dimensional 
spin systems. This equivalence could potentially provide 
new perspectives by enabling the use of methods and re- 
sults from one area in the other. A similar idea of baker's 
map construction has been used for establishing a con- 
nection between low-dimensional dynamical systems and 
universal computation 

Consider the following stochastic baker's map 



f{x,y)^{2x-\2x\,y/2 + £,{x,y)/2) 



(1) 



where ^(a;, y) is a stochastic variable that can take the 
values or 1, and where \x\ denotes the integer part 
of X. If ^ = [22;], we get the standard non-stochastic 
baker's map. In this paper we will consider a general 
case where ^ is characterised by a probability distribution 
that depends on the position (x, y) in state space. 

The dynamics of the baker's map is conveniently de- 
scribed using the dyadic representation (i.e., binary ex- 
pansion) for positions in the unit interval, using the no- 



tation X — {.xoXiX2-- ) and y — (.?;o2/iy2---)i so that x = 
X^i^o with Xi G {0, 1}, and similarly for y. The 

original non-stochastic map shifts the x sequence to the 
left removing and the y sequence to the right putting 
the symbol in the O'th position, i.e., x' = {.xiX2X^...) 
and y' = {.x^y^yi...). The two semi-infinite sequences 
can be put together {...y2yiyQ.xoXiX2---)- Using this rep- 
resentation the map is simply an operation shifting the 
combined sequence one step to the left. 

When we generalize the symbolic dynamics description 
to the stochastic baker's map, Eq. ([1]), an iteration again 
means that the xq symbol is removed but that the value 
of the stochastic variable ^ enters in the O'th position of 

y, 

f{---y2yiyo-XoXiX2...) = {...yiyo£,.xiX2X3...) (2) 

In general we will consider the case when the probability 
distribution for ^ depends on position [x, y), i.e., the full 
X and y sequences, and we denote the probability for ^ 
being by qo{x,y). 

In this paper, we are primarily interested in the case 
when the stationary measure of the map in the unit 
square is symmetric under exchange of x and y, since 
such a symmetry will be required when we make the con- 
nection to spin systems below. One way to achieve this 
is to consider a map y) that randomly chooses, with 
equal probabilities, either the function f{x,y) in Eq. ([1]) 
or the function f(jj, x). 



( f{x, y) with probability 1/2 
\/(j/,a;) with probability 1/2 



(3) 



The stochastic baker's map (j) is characterised by an in- 
variant measure over state space (x,y) together with 
an entropy rate s^^ of the map. The entropy rate can be 
derived using the symbolic dynamics approach p?]. The 



2 



following states illustrate three iterations of the map, 

■■■VnVn-l---VlV(3 ■ XiiXi...Xn^iXn... 

...?;„_iy„_2...yo6 
...y„?/„_i...?/iyo^6 

-■■y7i+iyn---y2yi^^3^2Xl...Xn^2Xn~l--- 



(4) 



where ^i, ^2, and ^3, are binary symbols introduced by 
the stochastic variable ^. Note that the probabilites for 
these may depend on the complete x and y sequences. 
Assume that we observe the system at a finite level of 
resolution with a binary precision n+1 as above. Then, if 
we have an infinite sequence of iterations, we will almost 
always know, by looking at the history, what is shifted 
in from lower levels of resolution. The only uncertainty 
comes from the stochastic component ^ and from which 
direction the sequence is shifted. The latter term just 
contributes with a constant, log 2. Assuming that the 
system is ergodic, which is supported by the calculations 
at the end of the paper for the stochastic characteristics 
we will use, the ergodicity theorem implies that a spatial 
average using can be used instead of a temporal aver- 
age. Then the average rate of entropy production can be 
written 



fib{x, y) a {qo {x, y)) + log 2, 



(t{p) 



(5) 



-plogp 



where cr(p) is the entropy function: 
(l-p)log(l-p). 

Next, we present the representation of the two- 
dimensional spin system that has a structure identical 
to the one we have used for the stochastic baker's map in 
Eqs. ((21 E]). The equilibrium state in a two-dimensional 
spin system with nearest neighbour interactions can be 
characterized by probabilities of spin configurations of 
the form in Fig. [T] 0, i, i, [13, [H 111 . 
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FIG. 1: The spin block configurations needed to capture the 
statistics that determines the equilibrium properties of a two- 
dimensional spin system. 



Denote the configuration of size 2n excluding the spin 
(, by Bniyn-iyn-2-yo-XoXi...x„-i) or shorter Bn{x,y). 
Spin variables are and 1. The conditional probability 
for spin ^ being given the configuration is denoted by 
Pn{x,y), and it is determined by the translation invari- 
ant measure fi that characterises the equilibrium system. 
In order to capture the full characteristics of the spin 
system, one needs to consider the limit n — > cxd. 

In this representation the relation between the stochas- 
tic baker's map and the conditional probability charac- 
terising the spin system is clear. By applying the con- 
ditional probability pn{x,y) (in the limit of n — > 00) we 



shift out symbol xq from the x sequence and we shift in 
symbol ^ to the y sequence. In this way we get a new 
spin configuration in the same way as we get a new po- 
sition in the baker's map. This leads to the main result 
of this paper: If we let the probability distributions for 
^(x, y) be identical to the conditional probabilities in the 
spin system. 



qo{x,y) = lim Pn{x,y), 



(6) 



we get a dynamics of the stochastic baker's map with 
an invariant measure /ifc that is equal to the translation 
invariant measure of the spin system, — fj,. 

The measure /i determines the statistical mechanics 
properties of the spin system. For finite block size n, let 
fin{x, y) denote the corresponding measure. Spin system 
symmetry requires ^n{x, y) = x), and it is reflected 

in the symmetric construction of the stochastic baker's 
map in Eq. ([3|). The entropy s (in terms of Boltzmann's 
constant fc^) of the spin system is given by 



s = lim y 

x.,ye{0.1}" 



Mn(a;,y)cr(p„(a;,y)), 



(7) 



Because of the construction of in the stochastic baker's 
map, by designing it from the characteristics of the spin 
system, Eq. the spin system entropy is up to a con- 
stant identical to the entropy rate of the baker's map, 
Eq. (O, s = Sf_i^ ~ log 2. The difi'erence between succes- 
sive improvements of the entropy estimate by increasing 
block size n can be interpreted as correlation information 
in the system 

Let us consider the Ising model with nearest neighbour 
interactions with interaction constant J = I (in terms of 
ks)- Then the energy can be written 

u = 2/11 (0, 0)(1 -2pi(0, 0)) - 2piil, 1)(1 - 2pi(l, 1)) (8) 

where /j,i(0, 0) denotes the probability for a Bi block to 
be (0.0), and similarly for the other term. The equi- 
librium state is characterized by a minimum in the free 
energy g (in terms of fc^) with respect to variations in 
the measure /x, 



u — Ts. 



(9) 



We now discuss a procedure that gives an approxima- 
tion of the invariant measure of the stochastic baker's 
map, and consequently it also determines the translation 
invariant measure fi of the spin system. 

If we choose qo{x,y) based on statistics from Monte 
Carlo simulations and run the stochastic baker's map (j>, 
does the system converge to a stationary measure that re- 
produces the physical properties of the spin system? One 
way to investigate this numerically, is to make approxi- 
mations by choosing a certain resolution for positions in 
the unit square in the baker's map as well as a certain 
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level of resolution for the dependence of qo{x, y) on posi- 
tion {x,y). For the positions in the dynamics, we divide 
the unit square into 2™ equal squares, which means that 
we use a binary precision of m. This is the resolution 
we will use to estimate the invariant measure of the 
baker's map dynamics. Since the conditional probabil- 
ity qo{x,y) is based on the statistics over the blocks of 
2n + l spins in Fig. [H the binary precision in the spatial 
dependence of qo{x,y) is n. We require that n < m. 

When we consider the baker's map approximated to a 
certain spatial resolution m it turns into a Markov pro- 
cess. The transitions can be represented by a finite reso- 
lution version hm of /, cf. Eq. ([2]), 

hmiym~l---yiyO-XQXi...X„i^2Xm-l) = ^^^^ 
= {y,n~2---yo^-XlX2-.-Xm-lO 

Here a new stochastic variable C enters at the finest level 
of resolution in the expanding dimension since a new 
binary symbol is shifted in from indistinguishable posi- 
tions. In order to shift in symbols in a consistent way, we 
use the available information that was used to form the 
conditional probability go- This gives us the conditional 
probability zo for C to be 0, given {x 

m — n- ■■•^m — l 

). We 

include this conditional probability in hm, and then we 
get a revised stochastic baker's map rjm by combining 
with its mirror image, as in Eq. ([3]). 

This means that we can calculate the invariant mea- 
sure at resolution m by considering a " master equa- 
tion" for the density function Pm{x, y; t) of the stochastic 
baker's map at the resolution m, 

Prn{x,y;t + 1) = ^ pm{x',y';t)P{x',y' ^ x,y) 

a;',a'e{0,l}" 

(11) 

where the transition probabilities P are based on the map 
r]m from Eq. ifTO]) and the discussion above. Using the 
double sequence representation, Eq. pT|) can be written 

Pm{yn--yi-xi...Xn;t + 1) = 

= ^ ^ Pni{y7i+l--y2-X0...Xn-i;t)x 

X qy^{yn+i...y2.XQ...Xn-l)Zx„{xo...Xn-l) + 
+ ^ ^ Pm{yn-l---yO-X2---Xn+i;t)x 

X q^^{xn+i.:X2-yo:-yn-i)zy„iyo--y7i-i)- 

(12) 

Here we use the notation gi = 1 — go and zi — 1 — zq. 

In other words, the map rjm defines a Markov process 
with a transition matrix P of size 2^™ x 2^™, and the in- 
variant measure fi^^ is then the dominating eigenvector 
of the P matrix. 

We illustrate the construction of stochastic baker's 
maps by using the procedure above for two temperatures 
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FIG. 2: The invariant measure ni, of the stochastic baker's 
map derived at resolution m = 8 and at temperature T = 4, 
calculated as the largest eigenvector of Eq. (|12p . The positions 
{x, y) in state space are numbered according to the binary 
expansion of the coordinates, i.e., from to 255. 

in the Ising model. The conditional probabilities for ^ 
and C, i.e., q and z in Eq. p2p . are in the first case, 
for temperature T = 4, based on statistics from Monte 
Carlo simulations on a 200 x 200 lattice, in total 3.1 x 10^ 
block configurations (as in Fig. [T|) with n = 3. The free 
energy is calculated by averaging the energy and esti- 
mating the entropy as in Eqs. (|7l8p . and the result is 
gMC ~ —3.03643 (same as the exact value at this pre- 
cision 13]). Then we solve the invariant measure 
for the corresponding baker's map, using Eq. p2p with 
a resolution m = 8. The invariant measure shown in 
Fig. [21 and its projection to the x axis in Fig. [31 both 
exhibit a fractal stucture as one may expect. We can 
now apply the invariant measure pb to Eqs. (|7l8p to cal- 
culate the free energy of the baker's map. The result is 
gb w —3.03634, which deviates less than 10~^ from the 
Monte Carlo value. 

0.025 p 
0.020 



0.015 
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FIG. 3: The invariant measure of Fig. [2] projected onto x. 
Positions x axe numbered according to the binary expansion 
of the coordinates. 

In order to get a qualitative picture on how the dif- 
ferent levels of resolution influence the construction of 
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the map, we design a series of maps based on different 
values for m and n. We choose a lower temperature, 
T = 2.2 just below the critical temperature Tc = 2.269, 
since correlations are stronger here and that may show 
up in the dependence on block size n. The construction 
of the baker's maps is now based on statistics from MC 
simulations, in total 3.9 x 10'' spin blocks with n = 5. 

The result is summarized in Table U showing that at 
sufficiently large blocks {n > 4) we get baker's maps that 
reproduce the free energy value of the MC simulation. 
Note that the spatial resolution of the state space (given 
by to) does not influence the free energy of the map as 
much as the choice of dependence on block size n. The 
invariant measure for to = 8 and n = 4 is shown with 
the projection to the x axis in Fig. HI Since we are below 
the critical temperature the left-right symmetry is bro- 
ken. This exercise illustrates that the invariant measure 
of the stochastic baker's map, also at finite resolution, 
is sufficiently close to the measure of the spin system to 
reliably reproduce the free energy. 



TABLE I: Free energy for the stochastic baker's map at dif- 
ferent levels of resolution of state space m and block size n 
for the conditional probabilities. The rightmost column shows 
the Monte Carlo estimates. The exact value of the free energy 
with this precision is —2.0907 13]. 



approach to find the invariant baker's map measure that 
minimizes the free energy. Such a procedure would have 



n 


m — 1 m = 2 


m — 3 


m = 4 


m = 8 


MC 


1 


-2.0807 -2.0804 


-2.0805 


-2.0808 


-2.0815 


-2.1009 


2 


-2.0894 


-2.0890 


-2.0888 


-2.0887 


-2.0928 


3 




-2.0905 


-2.0904 


-2.0902 


-2.0913 


4 






-2.0907 


-2.0906 


-2.0909 


5 








-2.0906 


-2.0907 



We have demonstrated that there is a class of stochas- 
tic baker's transformations that reproduce the equilib- 
rium measure in two-dimensional spin systems with near- 
est neighbour interactions. The generalization to longer 
(but still finite interactions) is straight-forward. This 
means that for any two-dimensional spin system, there is 
a corresponding stochastic baker's map that contains the 
equilibrium properties of the spin system. The difference 
between our result and those obtained for coupled map 
lattice (CML) models, e.g., [1], is that, in our approach, 
spin configurations (of a certain shape) are represented 
by positions in the unit square. This leads to the result 
that the process of moving across a spin lattice captur- 
ing the average physical characteristics is equivalent to a 
certain baker's map dynamics. In the former work, spa- 
tial configurations in the spin system are represented by 
corresponding spatial configurations in the CML system, 
and that cannot be used to establish the connection be- 
tween statistical mechanics and dynamical systems that 
we find in the present study. 

An interesting question is whether the finite resolu- 
tion map of Eqs. (|10llip can be used in a variational 
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strong connections to variational methods 
The reformulation of the variational problem for the spin 
system into a dynamical systems problem may offer new 
perspectives based on dynamical systems theory. Work 
along these lines is in progress. 

The author thanks Martin Nilsson Jacobi for several 
constructive discussions, and both him and Kolbj0rn 
Tunstr0m for valuable comments on the manuscript. 
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FIG. 4: The invariant measure /^b at T = 2.2, projected onto 
X, at resolution level m — 8. The first column, representing 
a; = is cut off as it reaches 0.55. The temperature is be- 
low the critical level Tc = 2.269, as is seen from the broken 
left-right symmetry, i.e., the symmetry between x and 1 ~ x. 
(Positions are numbered as in Fig. |31) 
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